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The influence of surface defects on the critical properties 
of magnetic films is studied for Ising models with nearest- 
neighbour ferromagnetic couplings. The defects include one 
or two adjacent lines of additional atoms and a step on the sur- 
face. For the calculations, both density-matrix renormaliza- 
tion group and Monte Carlo techniques are used. By changing 
the local couplings at the defects and the film thickness, non- 
universal features as well as interesting crossover phenomena 
in the magnetic exponents are observed. 
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I. INTRODUCTION 

Critical phenomena of magnetic films are of-cui-rent in- 
terest, both experimentally and theoreticallytrEl. In the 
limiting cases of one layer and of infinitely many layers, 
one deals with two-dimensional magpfitsQ and with stan- 
dard bulk and surface magnetismE EHl, respectively. For 
systems consisting of a finite number of layers, interest- 
ing crossover phenomena between these limiting cases are 
expected. 

In this article, we shall consider critical properties of 
ferromagnetic films of Ising magnets with various imper- 
fections at the surface, motivated partly by possible ex- 
perimental realizations of magnetic thin filrrisjaiitlj stripes 
of magnetic adatoms and stepped surfaceallO'tHi, partly 
by genuine theoretical interest. Imperfections may be 
due to regular or irregular changes in the surface cou- 
plings or due to additional structures on the surface. A 
simple example of the first case is a ladder of modified 
couplings in an otherwise uniform two-dimensional sys- 
tem as introduced by Barievt^, see Fig. 1(a). We will 
study this briefiy, since it can serve as a testing ground. 
Our main interest, however, is in additional structures, 
as depicted in Figs. l(b)-(d). Thus we will investigate 
surfaces with magnetic adatoms in the form of 

• one additional straight line. Fig. 1(b), 

• two neighbouring lines. Fig. 1(c), 



for various local couplings at the defects and for films of 
varying thickness. 
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a straight step of unit height. Fig. 1(d), 
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FIG. 1. Geometry and interactions of Ising models with 
various surface imperfections: (a) ladder defect, (b) one ad- 
ditional line of spins, (c) pair of adjacent lines of magnetic 
adatoms, and (d) straight step of monoatomic height. Usu- 
ally, surface spins are denoted by open circles, and bulk spins 
by full circles. In case 1(b), the shaded spins allow for both 
interpretations. 

Previous related work on Ising models includes the 
study of the step magnetization at the ordinary tran- 
sition of rather thick filmsEj and jtJac study of magnetism 
in thin films with rough surfacesEj. 

In studying the influence of these imperfections espe- 
cially on the critical behaviour, we use the denpJ», ma- 
trix renormalization group technique (DMRG)t£HlIl, be- 
ing most suitable in the case-pf merely one layer, and the 
Monte Carlo (MC) methodE3, which allows to treat films 
of considerable thickness as well. 

The article is organized as follows. In the next section, 
we present our findings on single layers with defects, ap- 
plying DMRG. The MC results on single layers and on 
films with an additional line of magnetic adatoms and 
with a straight step on top of the surface are discussed 
in Section 3. A short summary concludes the article. 



II. ONE LAYER: DMRG 

The planar Ising model with line-like defects is a pe- 
culiar system, because it shows non-universal magnetic 
exponents. This is connected with the values v ^ \ 
and Xs = 1/2 of the exponents for the correlation length 
and the surface magnetization of the pure system, respec- 
tively. A one-dimensional, energy-like perturbation then 
is marginal and can change the critical behaviour contin- 
uously. For this reason, the system has been the topic of 
various studiea3, -wjith the focus most receritly on a con- 
formal treatmentEj and on random systemsEJ. While the 
simple chain and ladder defects considered by Bariev are 
solvable free-fermion problems, the other cases we study 
are not integrable and one has to use numerical methods. 
In the following we dicuss the quantity of direct physi- 
cal interest, the local magnetization at or near the defect 
lines. 

To obtain it, we used the transfer matrix running along 



the direction of the defect, see Figs. l(a)-(c), and detet- 
mined its maximal eigenvector via the DMRG methodEll. 
In this way one is treating an infinitely long strip of 
width M with the defect located in the middle. Only 
the infinite-system algorithm was used, in which one en- 
larges the system step by step and always chooses an 
optimal reduced basis via the density matrix. This is 
very convenient, since one can insert different defects af- 
ter the system has reached the desired size. No further 
sweeps to optimize the state were made, since tests on 
the ladder defect gave very good results without them. 
Most calculations were done with 64 kept states and a 
truncation error around 10~^^. The local magnetization 
m(i) was determined from the spin correlation function 
C(i) — (cTiCTi) for free boundary conditions, or directly 
as {pi) for fixed boundary spins. The width was always 
much larger than the correlation length and varied be- 
tween M — 100 and M = 5000 for the temperature 
range studied (0.001 < t < 0.1, where i = 1 - T/T^ 
is the reduced temperature). The (absolute) error in m, 
determined by comparing with anaytical results was at 
most 10"''' for a system at t = 0.001, cut in the middle by 
a ladder defect. For less severe modifications and larger 
values of t it was even smaller. 
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FIG. 2. Spin correlation function C(z) for a strip of width 
M = 150 with ladder defects (below the plateau) or one ad- 
ditional line of spins (above the plateau), as obtained from 
DMRG calculations at the reduced temperature t — 0.072. 
The defect strengths Ji/Js and Ja/Js are indicated. Upper 
part: total view, lower part: central region. 

In Fig. 2 we show the correlation function C{i) across 
the strip for ladder defects (Fig. 1(a)) and for an addi- 
tional line (Fig. 1 (b) ; in the DMRG study we considered 
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the case J„ = Jg). The upper part gives an overall pic- 
ture, while the lower one shows the defect region in more 
detail. For ladder defects the strength J; of the defect 
bonds was varied, whereas for an additional line it was 
the coupling Ja between the line spins and the substrate. 
Since C{i) factors for large distances, these curves also 
give the profile of the magnetization in the bulk. One can 
see how m increases or decreases near the defect, depend- 
ing on the sign of the perturbation (similar curves were 
obtained inEJ for a random system) . If one cuts the lad- 
der bonds by choosing J/ = 0, one obtains the boundary 
magnetization of the homogeneous model in the middle 
of the strip. On the upper side, the possible increase of 
m depends on the details of the defect. It is limited if 
one varies J a, because a line with infinite J a is equivalent 
to a chain defect in the plane with merely doubled bond 
strength. 

The temperature dependence of m is shown in Fig. 
3 for the spins in the plane situated below one or two 
additional lines. One can see how it is increased over 
the Onsager value by increasing the coupling J a- As ex- 
pected, the effect is even stronger for two additional lines. 
In this case, m has already twice the undisturbed value 
for the smallest shown t. Quantitatively, this enhance- 
ment is described by a decrease of the exponent /3/, the 
local critical exponent which describes the vanishing of 
the magnetization near the additional line of magnetic 
adatoms. 
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FIG. 3. Local magnetization m of the spins below one or 
two additional lines as a function of temperature, for three 
values of the coupling ratio Ja/Js- The lowest curve is the 
Onsager result for the perfect Ising model. 

To investigate this, we have analyzed the temperature 
behaviour of ni mX&iaas of an effective (critical) exponent 
/3e//, defined byH'tile 

Peffit) ^ \n{m{U)/m{U+{))/\u{U/U+i) (1) 

with t = {ti+ti+i)/2 (alternatively, one could choose t to 



be the geometric mean t = ^titi+i). As one approaches 
the critical point, t ^ 0, this quantity converges to the 
true local exponent (3i . It is also a very sensitive indicator 
for the numerical accuracy of a calculation. 
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FIG. 4. Effective exponent /3e// as function of the reduced 
temperature t for one additional line and three different cou- 
pling ratios Ji/ Js- Full: Spins located below the line, dotted: 
spins in the line. 

Some typical results are given in Fig. 4 for one ad- 
ditional line and four values of the ratio Ji/Js of the 
couplings in the line. For J; — 0, one is treating the 
plane with independent attached spins and the Onsager 
result /? = 0.125 is recovered with high accuracy. In the 
other cases, the exponents both for the spin in the line 
and the one below it are shown and one sees that the 
two curves have different slopes, but a common limit for 
t — > which can be determined very precisely. The val- 
ues for (3i found in this way are accurate to at least three 
digits. For the case Ja/Js ^ 1 which, as mentioned, is 
equivalent to a line defect in the plane, this was checked 
explicitely by comparing with the analytical result. In 
the figure, also a negative J; is shown, which leads to a 
reduction of m and an increase of (3i over the Onsager 
value . In this limiting value 0.142 is approached 

rapidly for Ji/Js < —1. This is the same effect as for a 
chain defect in the plane with strong antiferromagnetic 
couplingsQ. In that case, the exponent is increased up to 
the value 0.5 of the free surface. The sign of Ja, on the 
other hand, has no influence on the exponent. 

The results for /?/ are collected in Table 1 and in Fig. 5, 
where the exponent is plotted as a function of the varied 
couplings (keeping the other couplings fixed and eaual 
to Js). For comparison also the analytical resultsE^B, 
for simple chain and ladder defects are shown in Fig. 5. 
One notes that, for a single line and small modifications, 
it does not matter much whether one changes Ja or J; . A 
large Ji/ Js, however, has a much more pronounced effect 
than Ja/Js, since it corresponds to additional spins which 
are almost rigidly locked together. For the double line. 
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the exponent drops much faster, reaching 10~^ aheady 
around Ja/Js ~ 1- For more additional hnes, i.e. for a 
terrace on the surface as in Fig. 1(d), this effect would 
be even stronger. In this case, the magnetization would 
practically jump as in a first-order transition. 
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FIG. 5. Local exponent j3i as a function of the ratio of the 
coupling strength A, as defined in the figure for five different 
situations. For chain and ladder defects the analytical results 
are shown, otherwise the DMRG results are depicted. 



TABLE I. Numerical values for the local exponent Pi of an 
Ising plane with one and two additional lines of spins. 
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III. FILMS: MONTE CARLO SIMULATIONS 
A. One additional line of spins 

Extending the DMRG calculations on an Ising layer 
with one additional line of spins, we did Monte Carlo 
simulations on the corresponding Ising films, consisting 
of L layers with one line of magnetic adatoms on top of 
the surface, see Fig. 1(b). We set J/ = Js, with Ja — 
Jn = Js (variant A, treating the spins directly below the 
additional line as surface spins, as it was done in the 
DMRG study) or Ja ~ Jn = Jb (variant B, treating the 
spins directly below the additional line as bulk spins). - 
In the layers, periodic boundary conditions are used. 

Let Sij^k = il be the spin on site in the kth 

layer. Taking layers of Af x spins, the spins in the 
additional line on top of the surface through the center 
are located at {i = {M + l)/2, j, k = 0), M odd, with j 
running from 1 to N. We computed, among others, the 
line magnetization m{i, fc; L) defined by 



'm{i, fc; L) 



(2) 



The magnetization of the line on top of the surface, m;, 
is given by mi{L)= m{{M + l)/2,fc = 0;L). 

In the simulations, the film thickness L ranged from 
1 to 40, with layer sizes being sufficiently large to cir- 
cumvent finite-size effects (up to 161 x 320). To speed 
up computations, the single-cluster-flip algorithm was 
implemented. We studied the cases (i) Jg = Jb as well 
as (ii) Js — 2Jb (variants A and B), which lead to the 
two characteristic scenarios of surface critical phenomena 
for L,M,N oo (semi-infinite case). In the first case, 
bulk and surface spins order simultaneously at tempera- 
ture Tc (ordinary transition) , while in the second one the 
surface spi|as. order at a higher temperature, (surface 
transition)LlJ3. 

(i) At the ordinary transition of the semi-infinite Ising 
model, L — + oo, the magnetization deep in the bulk 
vanishes lik£,r^ oc , with t = \T — Tc\/Tc, where 



(3 = 0.31...E3o. At the pesfeG.t, flat surface, one finds 
m <x tP\ with (3i k. O.SOEaM The vanishing of the 
magnetization in the additional line of spins on top of 
the surface is expected . to b e governed by /3i as well, 
i.e. /3/(L oo) = /3iBN. On the other hand, for 
a single perfect layer, L = 1, it is well known that 
TO oc 

ifed^ ^2<i = 1/8. Adding a row of spins, we ob- 
tain, from the DMRG calculations, to; oc 
l3i{L = 1) « 0.084, see Table 1. 

To monitor the influence of the layer thickness L at the 
ordinary transition, we computed magnetization profiles 
m{i, k] L), the critical temperature Tc{L), and the critical 
exponent Pi{L). The dependence of the transition tem- 
perature|-oii. the thickness L has been studied before for 
flat fllmsEltJ, and it is, certainly, not affected by the pres- 
ence of the additional line. In Fig. 6, the magnetization in 
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the defect line, mi{L), is depicted as a function of temper- 
ature for L ranging from 1 to 5, illustrating the increase of 
the transition temperature with L. In the ordered phase, 
T < Tc, the line magnetization m(i, k; L) is, in each layer, 
maximal for the center line, m{{M + l)/2, k] L), see also 
Fig. 2. The maximum is most pronounced at fc = 1 (we 
shall denote the magnetization in that line beneath the 
additional row of spins by mib — m{{M + 1)/2,1;L)), 
due to the increased coordination number, compared to 
the other surface lines. The magnetization in the addi- 
tional line, m;, is suppressed compared to wjb, because 
of missing neighbouring spins. 



m 




FIG. 6. Simulated magnetization in the additional line 
mi{L) for Ising films with L layers and one additional line of 
spins on the surface, choosing Js — Jt, for L ranging from 
1 to 5. Each layer consists of 81 x 80 spins. The critical 
temperatures Tc{L) are marked by arrows. 

Various crossover effects show up in the effective ex- 
ponent l3eff{i,k;L), defined by m{i,k;L) oc 
corresponding to the slope in a standard log-log-gl 
the temperature dependence of the magnetizationD'Ej 
see eq. (1). On approach to Tc, one expects to observe 
the limiting cases (3ef f{i, k; L) — > /3 for k and L large, 
— > /3i for k small and L large, — > f3i{L = 1) for L = 1 
and i = {M + l)/2, fc = or 1, and f32d for L = 1 
and sufficiently far away from the additional line in the 
centre. 

The crossover behaviour is illustrated in Fig. 7, show- 
ing /5e//((M + l)/2,fc;L) for m,(L), /c = 0, and m,;,(L), 
k — 1, with the film thickness ranging from L = 1 to 
10. (3eff{{M + 1)/2,0;L) decreases monotonically, ex- 
cept for L = 1, over a wide range of temperatures on 
lowering t, but with the effective exponent, at fixed i, 
increasing clearly with the film thickness, as depicted in 
Fig. 7(a). The data seem to indicate that the asymptotic 
critical exponent Pi{L), as t — > 0, of the magnetization in 
the additional line of magnetic adatoms increases, how- 
ever, only weakly with L, being quite small, around 0.1, 



for L going up to 10 (the increase itself may be argued 
to reflect the diminishing role of the defect line on the 
two-dimensional critical fluctuations in thicker films; of 
course, P{L) is bounded by 1/8 for finite L). For the mag- 
netization beneath the additional line, mih, corrections 
to the asymptotics are rather large as well, see Fig. 7(b). 
Here the effective exponent /3ef f{{M + l)/2, 1; L) changes 
with temperature in a non-monotonic fashion, except for 
L — 1. In agreement with the observations for mi, the 
true critical exponent f3i{L) is rather small, around 0.1, 
increasing only weakly with L. The location of the maxi- 
mum in (3ef f{{M + l)/2, 1; L) indicates the temperature, 
at which one crosses over from the regime dominated by 
two-dimensional critical fluctuations, close to the phase 
transition, to the regime, further away from the critical 
point, where the fluctuations are (nearly) isotropic and 
three-dimensional. Thence, at the maximum the cor- 
responding correlation length is argued to be about the 
thickness of the film L. In the thermodynamic limit, 
L ^ oo, the maximum is believed to shift towards t — 0, 
with its height being /3/ = /3i « 0.80. 

Note that the strong corrections to scaling, as seen 
by the deviations of the effective exponents from their 
asymptotic values, may cause severe difficulties in ex- 
tracting the true critical exponents in simulations as well 
as in experiments. 

Similar crossover phenomena, now between P2d, Pi and 
(3, are expected to occur for the magnetization far away 
from the defect line, when varying the film thickness. 
This aspect, however, is of minor importance in the con- 
text of this study. 
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FIG. 7. Effective exponent Pef}{{M + l)/2, fc;L) of the 
magnetization (a) in the additional line of spins (k = Q), mi, 
and (b) in the line beneath (k = 1), mib, for the Ising model 
with equal couplings, as obtained from the MC simulations. 
The sizes of the layers are up to 161 x 320 spins. Only results 
which are essentially free of finite-size effects are shown. 

(ii) At the surface transition of flat Ising films, the 
surface magnetization vanishes, on approach to Tj, like 
TO oc t^^"*, independent of L. The critical exponent for the 
magnetization at the additional line of magnetic atoms 
on top of the surface, with TO/(;5) oc i'^', depends on 
the local couplings at that line, as seen from our DMRG 
results for L = 1. Indeed, the situation is similar to that 
of the edge magnetization at the surface traiisition, the 
edge corresponding to an extended defect linecSB, where 
non-universality holds as well. 

For Js = 2Jfc and Ja ~ Jn = Js, variant A, one ob- 
tains for a single layer, from the DMRG method, f3i{L — 
1) « 0.084, i.e. the value is below that of the perfect two- 
dimensional Ising model because of the increase in mi due 
to the additional line of spins, see Fig. 2. The value in- 
creases weakly with layer thickness, becoming in the limit 
of the semi-infinite system — oo) = 0.091 ± 0.002, 
as inferred from MC data for films with thickness L up to 
40, and reasonable extrapolations. Because the critical 
fluctuations in a film of finite thickness are ultimatively 
of two-dimensional nature, one expects a non-universal 
critical behaviour at the defect line with /?/ depending on 
L. Actually, the slight increase of the critical exponent 
with L reflects the impact of the bulk spins, which now 
tend to lower the magnetization in the defect line. 

For Js — 2J(, and Ja = Jn = Jb, variant B, both for 
single layers, L — 1, and films, the magnetization profile 
TO(i, fc; L) close to Ts is non-monotonic exhibiting a min- 
imum at the center line i = [M + l)/2, see Fig. 8. This 
minimum is due to the reduction of the couplings J„ at 
the defect below the value Jg elsewhere in the surface. As 



one goes deeper into the bulk, the magnetization profile 
smoothens, which can be readily understood. 

The critical exponent describing the vanishing of mi 
(or mil)) depends rather weakly on the thickness L of the 
film. For L = 1, we estimate from the MC data (3i{L — 
1) = 0.38 ± 0.01, i.e. a value above the Onsager value of 
the perfect two-dimensional Ising model resulting, from 
the decrease of the magnetization at the defect lineEl. The 
effective exponent decreases on approach to criticality, 
t — > 0, when considering mih, while it increases when 
considering m; , allowing to estimate (3i accurately. -From 
data for fairly thick films, L up to 40, we estimate /3j 
of the semi-infinite system to be l3i{L = oo) = 0.34 ± 
0.02. The slight change of l3i{L) with L for films of finite 
thickness is, again, believed to be due to the correlations 
of the spins at the defect line with the bulk spins, which 
affect (3i in such a way that it is non-universal when the 
critical fluctuations are of two dimensional character. 



m(i,k) 




FIG. 8. Magnetization profile m(i, k; L) at different 
depths k in films with L = AO layers and one additional line of 
spins on the surface. The simulation was done for 81 x 80 spins 
in each layer, couplings Js = 2Jb, Ja = Jn = Jb and temper- 
ature T — O.OBTc. The magnetization in the additional line 
itseff is mi ^ 0.218. 



B. Step 

Finally, we briefly report our findings for the critical 
properties of the step magnetization. A straight step 
is introduced (actually two steps, to allow for periodic 
boundary conditions) by adding half a layer of. magnetic 
adatoms to the surface of the magnetic filmtfl, see Fig. 
1(d). We discrimate two couphngs, Jg if both neighbour- 
ing spins are surface spins, and Jb otherwise. We consider 
the line magnetization of the spins at the step edge, mse, 
vanishing on approach to the transition as TOje oc t^"" . 

For Js = Jh, i.e. at the ordinary transition, one ob- 
tains Pse ~ 0.80 in semi-infinite Ising models, L oo, 
i.e. the same value as for the critical exponent of the 
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surface magnetization, as had been shown in a previous 
Monte Carlo study on thick Ising films 'with a stepllj, in 
agreement with analytical consideration^^. However, in 
thin films, the critical behaviour is quite different. In the 
simulations, for a single layer L — 1 plus half a layer, we 
find a critical exponent close to 1 /2 (its concrete value de- 
pends rather sensitively on a very accurate determination 
of Tc), i.e. a value close to that of the surface critical ex- 
ponent f3i of the two-dimensional Ising model (note also 
its robustness against randomness in the couplingsEZI) . 
This observation can be understood in the following way. 
One is dealing with a composite system displaying, as 
the layer size goes to infinity, two distinct phase transi- 
tions, one at the critical temperature of the Ising plane, 
kBTc{L = l)/Js = 2.269..., and one at the critical tem- 
perature of the double layer, ksTciL = 2)/Js ~ 3.2. 
Related, composite Ising models have been investigated 
beforeE312J, showing that on approach to the upper crit- 
ical temperature, where half of the system is disordered, 
the critical behaviour of the magnetization at the inter- 
face (i.e., here, at the step) is governed by the surface 
critical exponent. The same scenario is expected to hold 
for finite films with Tc{L+l) > Tc{L). However, the tem- 
perature region where this behaviour can be observed, 
will become smaller and smaller as L increases. 

At the surface transition, the same considerations are 
believed to be valid. Indeed, in the case Jg = 2J;,, we 
found f3se to be quite close to f /2 for a single layer, L — 1, 
plus half a layer. For trivial reasons, l3se{L = 1) = 1/2 
holds for Js ^ Jf,, when the bottom layer and the ex- 
tra half layer decouple with the step edge being the sur- 
face of a two-dimensional Ising model. In the thermo- 
dynamic limit, where Tc{L -I- 1) = Tc{L), so that the 
above decoupling considerations do not apply, we esti- 
mated from MC data for films with up to 40 layers, a 
value of Pse — 0.33 ±0.02. Presumably, in that limit, /3se 
is non-universal at the surface transition, depending on 
the ratio Js/ Jb- 

IV. SUMMARY 

Using density-matrix renormalization group and 
Monte Carlo techniques, we studied critical properties 
of magnetic Ising films with various surface defects. 

In particular, the effect of the local couplings at one or 
two additional lines of magnetic adatoms on the surface 
as well as at straight steps of monoatomic height has 
been investigated, especially in the limiting cases of films 
consisting of merely one layer and rather thick films. 

In the case of a single layer, L = 1, with additional 
lines of magnetic adatoms, the critical exponent of the 
magnetization at the surface defect is non-universal. The 
dependence of its value on the local couplings, as com- 
pared to that of the perfect two-dimensional situation, 
follows the trends observed for the exactly soluble two- 
dimensional Ising model with ladder and chain like bond- 



defects. The value may be lower or larger than in the 
perfect situation, 1/8, corresponding to an increase or 
decrease in the magnetization at the defect line. Adding 
half a layer of spins, one recovers, at the step, the sur- 
face critical exponent, 1/2, of the two-dimensional Ising 
model. 

In the limit L — > oo, varying the strength of the surface 
couplings may lead either to a surface or an ordinary 
phase transition. The change of the critical exponent of 
the magnetization at the defect has been found to depend 
only fairly weakly, for both types of transition, on the film 
thickness L in the case of one additional line of spins. 
At steps, the critical exponent is argued to be 1/2, for 
films of finite thickness and both kinds of transition, in 
agreement with the simulations. 

In the paper, we have not only presented the results 
for the exponents, but also shown various magnetization 
curves directly, so as to give an impression of the size of 
the effects. This is meant to encourage further experi- 
mental work on such surface structures and their mag- 
netic properties. 
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